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We prove that in three space dimensions a nonlinear wave equation utt — Ait = u v with p > 7 being 
an odd integer has a countable family of regular spherically symmetric self-similar solutions. 

1. Introduction 

An important feature of many nonlinear wave equations is that their solutions, corresponding to smooth initial 
data, may form singularities after a finite time. Such a phenomenon, usually referred to as blowup, has been a subject 
of intensive studies beginning with the pioneering works by Keller [l[ , John Q and Glassey 0] (we refer an interested 
reader to the excellent online review of this subject with the complete bibliography Q). 

In this paper we consider the semilinear wave equation with a power nonlinearity 

$ tt - A$ - $p = 0, $ = $(t, x), xeR 3 , (1) 

where p > 7 is an odd integer. The sign of the nonlinear term corresponds to focusing which means that it tends to 
magnify the amplitude of the wave. If $ is small this term is negligible and the evolution is essentially linear leading 
to dispersion. However, if <!> is large the dispersive effect of the laplacian may be overcome by the focusing effect of 
the nonlinearity and a singularity can form. Actually, neglecting the laplacian altogether and solving the ordinary 
differential equation <& tt = $ p one gets the exact, homogeneous in space, solution 



2(p+l) 



(P-1) 2 



a= T>0, (2) 

p-1 



which blows up as t — > T. By the finite speed of propagation one can truncate this solution in space to get a solution 
with compactly supported initial data which blows up in finite time. There is theoretical [j| and numerical Q evidence 
that the solution <f>o determines the leading order asymptotics of blowup for generic large initial data. However, for 
specially prepared initial data, in particular for data fine-tuned to the threshold for blowup, singularities may have a 
different form which is given by self-similar solutions of equation ([I]) . Such solutions were found numerically in [|| . 
The aim of this paper is to give a rigorous proof of their existence and discuss their properties. 
By definition, self-similar solutions are invariant under rescaling 

$(t, x) - * A (t, x) = \- a $(t/\, x/X), (3) 

hence in the spherically symmetric case they have the form 

<P(t,r) = (T-ty a u(p), P=j^~ t , (4) 

where T is a positive constant. Note that each self-similar solution, if it is regular for t < T, provides an explicit 
example of regular initial data developing a singularity in a finite time. Substituting the ansatz j4|) into equation ([1]) 
we obtain the ordinary differential equation for the similarity profile u(p) 



(1 - p 2 )u" + Q - (2 + 2a)pj u' - a{a + l)u + u p = . 



(5) 



It is easy to see that this equation has the constant solution uo(p) = bo which, of course, corresponds to the homoge- 
neous solution <&o of equation (Q]). As mentioned above, numerical results indicate that, besides uq, there exist also 
nontrivial regular solutions of equation ([5]). We remark that equation ([5]) has been studied by Kavian and Weissler 
who made many interesting observations about the behavior of solutions but unfortunately they imposed very 
restrictive fall-off conditions at infinity which excluded nontrivial solutions. However, as long as one wants to have 
an example of blowup, due to the finite speed of propagation, only the behavior of solutions inside the past light cone 



2 



of the blowup point (i = T, r = 0) is relevant, which corresponds to the interval < p < 1. As we shall see below 
regular self-similar solutions are parametrized by their values at the endpoints c = u(0) and b = u(l). Due to the 
singular nature of equation (JS]) at p = and p = 1, the generic solution obeying the regularity condition at p = 
becomes singular at p = 1 and vice versa, hence globally regular solutions can exist only for discrete values of the 
parameters b and c. The goal of this paper is to show that there exists an infinite sequence of pairs (b n ,c n ) which 
give rise to globally regular solutions. To this end, in section 2 we first prove existence of local solutions near the 
endpoints. Then, in section 3, we show, using a shooting argument, that for a discrete set of values of the parameters 
these local solutions match smoothly at a midpoint. In section 4 we show that the solutions constructed in section 3 
can be extended beyond the past light cone to infinity. Finally, in section 5 we derive some scaling properties of the 
shooting parameters. In the appendix we present an alternative rigorous proof of existence. 



2. Local existence 



In the first step we will analyze the behavior of solutions near the boundary points p — and p = 1 and prove the 
local existence of one-parameter families of regular solutions there. 

Near p = we introduce a new variable v = u' and rewrite equation (|5|) as the first order system 

pv! = pv, (6a) 

PV = -2- + 7—^2 {j— 1 P v + uP ) ■ ( 6b ) 

This system has the form required by Proposition 1 of Q guaranteeing the existence of a one-parameter family of 
regular solutions near p = such that 

<P) =c+l (i^^c \A p* + 0(p% (7) 



3\(p-l) r 2 

Similarly, near p = 1 we put s = 1 — p, u(s) = u(p), v(s) — u'(s) to get 



su — sv, 



1 (2 2(p+l)(l-s)\ 2(p+l) 



(2-sHVl-s p-l J (p-1) 



(8b) 



In order to bring this into the form required by Proposition 1 of 3 we introduce v = v — 2(^-1) u ^~£~ uP an< ^ rewrite 
the system (8) as 

2 

v + sf(u,v,s), (9b) 



p-l 



where the function f(u,v,s) is analytic in u, v and s near s — 0. We conclude that there is a one-parameter family 
of regular solutions near p—1 such that 



(10) 



The singular solutions behave like u ~ (1 — p)^ 1 near p = 1 leading to a singularity of u 1 . For later use it is convenient 
to introduce the variables a = (1 — p) 2 /^ -1 -*, ip( a ') = u (p) an d 0(a) — au' . The system (8) may then be written in 
the form 

V/ = -^-a^e, (11a) 

9' = g(i>,9,a), (lib) 

where the function g(ip, 9, cr) is analytic in ip, 9 and a near a = 0. This implies that there is a two-parameter family of 
regular solutions parametrized by the values b = ip(0) = u(l) and d = 9(0). These solutions are analytic as functions 
of 6, d and a in a neighborhood of a — 0. Thus, the general solution is in fact regular at p — 1 as a function of 
a = (1 — p) 2 /(p _1 ) with a convergent Taylor series at a = 0. 
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3. Global behavior and proof of existence 



The main result of this paper is: 

Theorem 1. For any odd p > 7 and any non-negative integer n, equation ([5]) on the interval < p < 1 has an analytic 
solution u n (p) which satisfies the boundary conditions ([7]) and (jlOp . and for which the function w n = u n /Uoo — 1; with 
Moo defined by (fl8|) . has exactly n + 1 zeros. 



Remark: In this section we give a "physicist" proof (no epsilon-delta stuff) based on scaling arguments. We believe 
that sacrificing slightly the mathematical rigour we gain better understandability of the basic mechanism which is 
responsible for the existence and structure of solutions. An alternative rigorous proof is presented in the appendix. 

Proof: First, we need to show that solutions which are regular at p = or p = 1 remain bounded on the whole interval 
< p < 1. This is conveniently proven using a suitable Lyapunov function [9j defined as 

ff-fl n*) U '\ UP+1 {P+1)U2 (12) 

obeying 

H{p)>-JL-(*<P+Q)* and fl<=(£±lp-Hy><o. (13) 



p — 1 \ (p — l) 2 / Vp — 1 

Hence, for solutions regular at p = we get 



C P+1 Cjn _i_ i~) c 2 

F(p)< J ff(0) = ^— -f^J- (14) 
p + 1 (p - l) z 



which implies 



\/r^7K(p)i < C ^, (15) 

and 

\u(p)\<c for c> (2(p+l)/(p-l) 2 ) 1 /( p - 1 ). (16) 

Furthermore, H{p) and thus w(p) has a finite limit at p = 1. 

On the other hand, for solutions regular at p = 1 we may integrate the function i/ to the left. It is easy to see that 

- (H + ly'H' < 2/p, (17) 

hence H stays finite on < p < 1 and with it u and u' . We emphasize that this argument does not exclude solutions 
which are regular at p — 1 and singular at p — 0. Actually, such an explicit singular solution 

i 

(2 3)\ p_1 
jp~~if) ' ^ 

— 2 

will play an important role in our analysis. Note that corresponds to the singular static solution $ = b^r p- 1 of 
equation |T]). 

Next, we consider the behavior of solutions regular at p — 0, i.e. solutions satisfying the initial condition (JT]), with 
large values of c = u(Q). To this order we rescale the variables 

u(p) = cU{x). (19) 



c (p-l)/2' 



From equation (O we get 

Pi + -f^ = ^( X ^ + 2 J^J^ + p^u). (20) 
ax x ax c v 1 V dec p — 1 dx (p — lj z / 



In the limit c — > oo we obtain the simple equation 
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on the interval < x < oo together with the boundary condition at x = 

U(x) = 1 -x 2 /6 + 0(x 4 ). (22) 

Note that the special solution it^ is invariant under the performed rescaling and thus Uoc(x) — frooX -2 '" 3-1 ) solves 
equation (I21|) . From our numerical analysis we know that U oscillates around Uoo suggesting the change of variables 

U(t)=x 2 ^-^U(x), r = ln(x). (23) 

Substituting (f23f into equation (f2Tj) we get the autonomous equation 

d 2 U , p-5dt) , - 2(p-3). 



dr 2 p — 1 dr (p — 1) 



t/ p - n C/ = 0. (24) 



Using the Lyapunov function 



ld?7 2 C/p +1 (p-3)J7 2 d/i p-5/dU\i 

h = H ; -r- With = 

2 dr p+1 (p - l) 2 dr p - 1 V dT / ' 



(25) 



we conclude that for r — > oo the solution t/ (r) tends to the fixed point = of equation (f!M|) and =p — > for 
r — » oo. Putting U = boc + y and ^ = zwe find 

g = z, (26a) 

^ p — 5 2(p — 3) ... . ,„„, . 

— = ~—:Z- Ky / y + N(y), (26b) 
dr p — 1 p — 1 

where iV (y) denotes non-linear terms in y. Neglecting these non-linear terms we get a linear system with the eigenvalues 
A = (— p + 5 ± i\/lp 2 — 22p — l)/2(p — 1). This implies that the solution of equation (f2T|) which is regular at x = 
has for x — > oo the asymptotic form 

1/ (a;) = ii^FT + Aoa; -5 ^^ sin(wlna; + 5 )J (27) 



with some constants Ao and 8q and = \Jlp 2 — 22p — l/2(p — 1). 

In order to prove the existence of regular solutions for sufficiently large c we choose some xq, such that for x > xo 
we can safely use the asymptotic expression (|27l) for the solution of equation (f2"Tj) with the boundary condition (|22[) . 
Furthermore, we choose Co so large that we can neglect the r.h.s. of equation (f2T))) on the interval < x < xq for 

c > Co- Integrating solutions regular at p = to po = xq/c^ and varying the parameter c = m(0) we get a smooth 
curve Co in the (u(po), u'(po)) plane. Using the asymptotic expression (|2T)l for the part of Co corresponding to c> c 
we get 

u(p) ps baop^^ (l + Aqc^ p^^T) sin(wln(c £ ^ i p) + ^ )J ■ (28) 

Keeping p = po <C 1 fixed and letting c — > oo (and with it a;) we find that Co spirals down (in the negative direction) 
to the limit point 

P(po) = (booPo 2/ip - 1 \ -^fcooPo (29) 

p—\ 

Having that it is easy to show existence of regular solutions for large c. Integrating the solutions regular at p = 1 
back to po we obtain a smooth curve C\ parametrized by b. For < b < bo this curve goes from the point (0, 0) to 
(&0j0) passing through P(po) for b = boo. Thus, the curve Co intersects the curve C± again and again as c goes to 
infinity (note that the curve C\ has positive rotation around P(po)). This is illustrated in Fig. 1 in the case p = 7. 
Obviously each such intersection yields a solution which is regular on the whole interval < p < 1. As c and thus the 
domain for x increases the corresponding solutions U(x) perform more and more oscillations around Uao^x). In fact, 
each further oscillation of U(x) about zero leads to another cycle in the (it, u') plane at p a . Thus, we arrive at the 
conclusion that for any sufficiently large n there exist regular solutions of equation {5| on the interval < p < 1 with 
n of such oscillations. 
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FIG. 1: The intersection of curves Co and Ci for po = 0.1 (p = 7). 

In the second step of our existence proof we shall show that we can decrease n and thus prove existence of regular 
solutions for any positive integer n. The constant solution uq with c = cq = bo has exactly one such intersection. 
Thus we must be able to decrease the number of intersections to one decreasing c from some large value c 3> 1 
corresponding to a regular solution with n»l intersections. Obviously the number of intersections can only change 
at p = 1, because, by the miiqucncss theorem for solutions of differential equations, a double zero of u(p, c) — Uoo(p) 
is not possible at a regular point < p < 1. On the other hand, no solution regular on the whole interval < p < 1 
can have b — u(l) = boa. Thus, a change of n can only happen for a value of c corresponding to a solution which 
is singular at p = 1. We showed in Section 2 that such solutions are parametrized by the numbers b = u(l) and 
d = V(l), where the latter is the finite coefficient of the singular mode. For solutions regular at p = the coefficient 
d depends smoothly on c. Suppose b(c) = boo for some c = c s . Then, we get a decrease of n exactly by one if 
either b — boo moves from positive to negative values as c decreases through c s and d(c s ) < 0, or b — boo moves from 
negative to positive values as c decreases through c s and d(c s ) > 0. This shows that as we decrease c from the value 
corresponding to a regular solution with k 3> 1 intersections to Cq we necessarily encounter regular solutions with any 
value of n between one and k. This concludes the proof of Theorem 1. 

Clearly, there could be more than one regular solution with the same number of intersections, because b — boo can 
have zeros of even order. However, our numerical analysis shows that there is exactly one regular solution for any n. 
The first few solutions generated numerically are shown in Fig. 2. 



4. Extension beyond the past light cone 



In this section we show that the solutions constructed above can be smoothly extended beyond p = 1 to infinity. 
Consider the function 



Q(P) = \(l P 2 )p 3 u' 2 + \p\\ - p 2 )uu> + 



3(5 -p) 



4(p - 1) (p - l) 2 



p 3 u 2 + 



1 



p+ 1 



p 3 uP +1 . 



(30) 



This function was introduced by Kavian and Weissler Q in their study of equation ((5]). By straightforward (but 
tedious) computation one gets 



2(3u/2 + pu') 2 + 3P 7 V + ^ U p +1 



2(p-l) 



P+l 



(31) 



hence Q(p) is monotone decreasing for p > 5. Since Q(0) = for analytic solutions, it follows that Q(p) < for p > 
and thus u(p) is positive for < p < 1. 

Remark. We note in passing that for p = 5 it follows from (f3"Tj) that the function Q(p) is the first integral which 
implies that uq is the only regular self-similar solution. 



6 




FIG. 2: The first five self-similar solutions for p — 7. To show oscillations around u«, we plot w n = u„/uoo — 1. 



It is convenient to rewrite equation J5]| in a self-adjoint form 



(/(l-p>')=/(l-P 2 ) 



2\a-l„ 



■PS" 



(32) 



Proposition 1. If it(p) is a regular solution of equation (|32|) defined for < p < 1 and w(0) > &o, then u(l) < bo. 

Proof: First, note that if u is analytic at p = 1 and it(l) = 6o, then u(p) = 6o, which contradicts the assumption 
that u(0) > 6o- Thus, the case u(l) = 6 is impossible. Next, suppose that u(l) > bo- If u(p) > bo on the whole 
interval < p < 1 then the r.h.s of equation (|32|) is negative and hence by integrating it we get that u'(l) = oo. 
Thus, u{p) < bo for some < p < 1 and therefore u must have a minimum value at some po, i.e., u'(po) = and 
< u(po) < fro- If so, between po and 1 there is a p\ where u(p\) = bo and u(p) < bo for po < p < pi- We now show 
that this is not possible. To this order we integrate equation (|32|) over the interval po < s < p < p±. Since u'(po) = 
we get 



p 2 (i- P 2 Tu>(p) 



s 2 (i- s 2 ) Q -y&[T 



)ds. 



S 2 (l-,s 2 ) c 



(33) 



We have u(s) < u(p) < b , hence u^' 1 - U?- 1 ) = u(b^ 2 + 6g" 3 M + ... + uP- 2 )(b Q -u)<(p- l)^" 1 ^ - ")■ Thus, 



uP-^ds < (p - ljftg-^fto - u(p)) / s 2 (l - .s 2 )"- 1 ^ = /(p)(6 - <p)) (34) 



for some bounded continuous function f(p). It follows from (153")) and (|3"4")) that 



< 



b -u(p) p 2 (l-p 2 ) 



2\a 



9{P), 



(35) 



where g(p) is a bounded continuous function on po < p < pi < 1. Integrating (|35[) from po to p, we get ln(6o — 
u(po)) — ln(&o — u {p)) < /X, 9i s )ds, hence ln(&o — u(p)) is bounded for p < pi. This contradicts the assumption that 
u(p\) = bo and concludes the proof. 

Next, we consider the behavior of solutions for p > 1. 

Proposition 2. If u(p) is a solution to equation (|32|) that is analytic at p = 1 and m(1) = 6 < b , then u(p) is defined 
for all p > 1 and < u(p) < bo- 
Proof: Integrating equation (|32[) from 1 to some p > 1 we get 



pV-i)V(p) 



1 (l-s 2 ) a u{b p - 1 ^uP-^da 



(36) 
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hence if u is bounded, so is u' and therefore the solution cannot become singular. Thus, it suffices to show that 
< u(p) < b. From equation (J32) we have that u'(p) = and < u(p) < b for p > 1 implies that u"(p) < 0, thus u is 
monotone decreasing. To show that u cannot become negative we introduce the function f(p) = pu'(p) + (a+l)u(p)/2. 
Evaluating /' at / = we obtain 



1 , „ 2p 2 
= -(1 -a 2 )u+ — H 
f=o 2 V ' p 2 



1 



-v? > 0, 



(37) 



which, together with /(l) = V /2a > 0, implies that f(p) > for p > 1. Since u'{p) < (as noted above), we conclude 
that u > 0, as claimed. 



5. Large n asymptotics 

Our numerical analysis shows that the shooting parameters b n and c n exhibit remarkable scaling properties in the 
large n limit. In this section we explain this phenomenon. 

In order to derive the large n asymptotic behavior of solutions which are regular at p = 1 we introduce a new 
variable w defined by u — Moo(l + w )- Using equations © and (|10[) we find that w satisfies the equation 

P 2 (l - P 2 )w" + {^^P 2p 3 ) w' - + «0((1 + v>7- 1 - 1) = 0, (38) 

with the boundary condition at p = 1 

w{p) = b _^. + _ W^){p - 1) + 0((p - I) 2 ). (39) 

For any < po -C 1 and & sufficiently close to feoo, the solution of equation (|38|) satisfying the initial condition (|39|) 
will stay arbitrarily close to w = on the interval po < p < 1, hence it can be approximated by w = fc ~ b °° wj,, where 
u>l satisfies the linearized equation 

p 2 (l - p 2 )< + (^~^P 2p 3 ) < - ^TJ-l = 0, (40) 

with the boundary condition wl(1) = l,w^(l) = (p — 3)/2. We could express tut in terms of a hypergeometric 
function but for our purposes it is sufficient to have its asymptotics for small values of p 



w L w p^=T)Aism(ujlnp + 5i), lu = y/7p 2 - 22p - 1/2 (p - 1) , (41) 
where Ai and 5i are constants. Thus, the solution regular at p = 1 with 6 near b^, has the following form for small p 

u(p) ~ Moo f 1 + & b °° p^T)Ai sin(wlnp + £i) J . (42) 



The key point is that the regions of validity of the approximations (|4"2")l and (|2"8")l overlap so we can match them. 
Matching the amplitudes of solutions (|28| and (142)) at po we get the condition 

A c^ =±A 1 b —^. (43) 

Ooo 

On the other hand, using (|28p and the fact that the phases of two adjacent solutions u n (p) and u n +i(p) differ by 7r 
we obtain 

hence 

« e<* , (45) 
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n 


Cn 


b„ 


Ac„ 


Ab n 


1 


2.054390385 


0.688698572 


2.8018 


0.4713 


2 


5.756037116 


0.820493408 


2.4090 


0.7428 


3 


13.86655615 


0.746908360 


2.5899 


0.5670 


4 


35.91343330 


0.796055093 


2.4577 


0.6752 


5 


88.26661166 


0.766263419 


2.5326 


0.6059 


6 


223.5507381 


0.785548198 


2.4823 


0.6493 


7 


554.9215495 


0.773546658 


2.5128 


0.6217 


8 


1394.439242 


0.781209584 


2.4930 


0.6391 


9 


3476.402010 


0.776393971 


2.5053 


0.6281 


10 


8709.676250 


0.779451184 


2.4974 


0.6351 


11 


21752.40861 


0.777522645 


2.5012 


0.6307 


12 


54434.14714 


0.778744142 


2.4993 


0.6334 


13 


136047.6759 


0.777972446 


2.5000 


0.6317 


14 


340293.1022 


0.778460765 


2.5008 


0.6328 


15 


850746.1358 


0.778152079 


2.5003 


0.6321 


oo 


oo 


0.778271716 


2.5005 


0.6324 



TABLE I: Shooting parameters of solutions u„ and their quotients Ac n = c n +i/c n and Ab n = (6„+i — 6oo)/(6oo — b n ) for p = 7. 
The last row for n — oo corresponds to the analytic results (45) and (46). 



which, using (|43f . implies that 

b n +l ~ boo - -»- B > f Af\\ 

— ~ — w e 2 (p-i)^ . (46) 

Oqg On 

The numerical verification of these scaling properties is shown in Table 1 in the case p = 7. 

Note added in proof. Peter Breitenlohner (personal communication) pointed out that a minor modification of our 
proof applies also to even values of p, hence the conclusion of Theorem 1 holds for all integers p > 6. 
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I. APPENDIX 



We present here an alternative version of the proof of Theorem 1. Throughout the appendix the range of p is < p < 1. 

We know from Section 2 and equation (16) that given any c there is a unique solution u(p, c) of equation (5) 
satisfying u(0, c) = c defined for all < p < 1. For p > we define v(p, c) = u(p, c)/uoo(p) and w(p, c) = v(p, c) — 1. 
We define also the function 

1 fv 2 vP +1 \ 

H v = V(l - p 2 )v' 2 - a(l - a) ( V - - ^— . (47) 



Using equation (5) we get 



2' 



K = ^pv' 2 , (48) 
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hence H v is decreasing along solutions. Furthermore, H v has a minimum equal to — (p — 3)/(p 2 — 1) for v' = and 
u= 1. 

Lemma 1 Let d = (p — l)/p. If c > 6 /d, then there is an p! = pi(c) < [boo/dc) 1 / 01 such that w(pi,c) = and 
< piw'(pi,c) < a. 

Proof: It follows from equation (32) that u(p, c) is a decreasing function of p for c > bo as long as u > &o- Integrating 
equation (32) and noting that < u(p, c) < c we get 

p 2 {l-p 2 ) a u'(p) = - [ s 2 (l-s 2 ) a - 1 [u(s) p ~a(l + a)u(s)]ds > - f s 2 (l-s 2 ) a - 1 u(s) p ds > - [ s 2 {l-s 2 ) a - 1 c p ds. 
Jo Jo Jo 

(49) 

Thus 

fP s 2 (l-s 2 ) 
(1-P 2 ) 



p 2 (l - P 2 Yu'(p) > - I ,) ' cJ>ds > - I s 2 cfds = - y -<3> , (50) 



hence 



u'{p)>- P <?. (51) 
3(1 - p 2 ) 



Integrating (|5ip from to a; where u(x, c) = dc gives 



- C (l-d)>|ln(l- 2 ; 2 ), so ^l-exp^--^ > . (52) 

Let y be a point where Uoo(y) = 6oo2/~ Q = dc, that is, y = [b^/dc) 1 ^ . We have 

boo\ 2/a _(b 00 y- 1 _( P y-^ip-s) e 



, / ■ 1,2 < - > ( 53 ) 

p — 1 / (p — l) 2 p 

hence x > y, and therefore u(y, c) > u(x, c) = dc = u oc (y), or w(y, c) > 0. Since w(p, c) < for small p, we conclude 
that there is a p such that w(p,c) — and, if p\ is the smallest such p, then w'(pi,c) > 0. Note that p\ < y. This 
conludes the proof of Lemma 1. 

Corollary. w(p) > —2a for p > p%. 

Proof: Using the fact that the function H v (p) is monotone decreasing we obtain 

H v (p) < H v ( Pl ) < ^ - a(l - a) Q - , (54) 

for p > pi- On the other hand, for p > pi 

H v (p) = ip 2 (l - p>' 2 - ail - «) ( j - ^) > -a(l - a) ( - - ^) . (55) 

Combining equations ([54")) and ([55]) we have 

a 2 / 1 1 \ (v 2 vP +1 \ 

f(v) = — - a{l-a)y-~- J-^jJ + «(1 -«)( T - R )>0. (56) 

for p > p\. Note that /'(«) > for v < 1. A straightforward calculation yields /(l — 2a) < if p > 7, hence v > 1 — 2a 
(or ui > —2a) for p > pi, as claimed. 

Let us define 



R{p,c) = y/w(p,c) 2 + w'(p,c) 2 , 0<p<l. (57) 

By uniqueness of solutions of ODEs, a solution starting at (u,u') — (c, 0) at p = cannot have {nip^)^' [p,c) = 
(woo(p), m oo(p))j which implies that i?(p, c) > 0. Thus, we may also define 

0(p, c) — arctan ( E^Lj£l j up to a multiple of 7r . (58) 
V w{p,c) J 
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Since the region {(p, c) 1 < p < l,c > 0} is simply connected we may unambiguously define a function 0(p, c) once 
we specify its value at any point in the domain. Since w(p, 0) = — l,w'(p, 0) = 0, we set 0(1/2,0) = ir. Note that 
lim p ^ + w (p, c) = — 1 an d lim p ^ + P w '(p, c ) = 0, hence lim p ^ + 0(p, c ) = tt- 

Lemma 2. Let p 2 <C 1. We have lim^oo 0(p2, c) = — oo. 

Proof: Consider the equation satisfied by the function 0(p) = 0(p, c): 



sin 2 9 + (P , ^vf^^ sin © cos 9 + + „ + ...+ cos 2 9 . 

(p-i)(i-p) (p-i) 2 (i-p 2 ) 



(59) 

We will show that the quantity in brackets in bounded from below by a positive constant r\ for p > pi . It follows that 



lim 6(p 2 ,c) = lim (e(p 1 (c),c)+ e'(p)dp) < J + Urn / 

c^oo c^oo y J pi(c) J 2 c-oo^ 



=-- v lim In = -oo , (60) 

p) 2 '«~oo Pl ( c ) V 7 



Pi( c ) 

and the lemma is proven. 

To get a bound on the quantity in brackets we regard it as a (formal) quadratic form in sin and cos0: [] = 
A sin 2 6 + B sin 6 cos + C cos 2 6, where 

(p-l)(l-/3 2 ) (p- 1) 2 (1 - p 2 ) 

Since p < p2, we may replace these coefficients by 

■4 = 1. e-^(i + . + ...+^). (62, 

We need to show that the discriminant A = B 2 — AAC is negative for v > 1 — 2a. We have 

Clearly, (p — 1) 2 A is a decreasing function of v. When w = 1 — 2a we get 

(p - 1) 2 A = (2p 2 - 8p + 6) -p 2 + 6p-5, (64) 



which is negative for p > 7, hence A < which concludes the proof of Lemma 2. 

Next, we consider solutions U(p,b) which start from initial value (10). As we showed in Section 3 these solutions 
arc defined for < p < 1, hence we may define 



R(p,b) = ^W(p,b) 2 + W'(p,b) 2 , 0<p<l, (65) 

where W(p,b) = U(p,b)/u OQ (p) — 1. By uniqueness of solutions of ODEs, a solution starting at p = 1 cannot have 
(U(p, b), U'(p, b) = (ii^fp), u'oo(p)), unless b = b^, which implies that R(p, b) > if b ^ b^. Thus, we may also define 

0(p, b) = arctan [ ^ ^— ' ^ ] up to a multiple of it . (66) 



VK(p,6) 

Since the region {(p, 6)|0 < p < 1,0 < b < boo} is simply connected we may unambiguously define a function 0(p, 6) 
once we specify its value at any point in the domain. We set 0(1, 0) = ir. 
We define a map 

<5>:R + = {c\c>Q}^R 2 + = {(x,y)eR 2 \y>0}, $(c) = (6(p ,c), R(po,c)) , (67) 

and a map 

,t, i C ( \ r2 lT , /.n / ( (Po, 6) - 2nn,R(p a ,b)) if < t < b^ (segment 1), 

I (7T - 2n7r, 1 — t) it t < (segment 2). 

Note that if ^(fr) = <&(c) for some b S (0, &oo) an d some c, then <I>(c) cannot lie on segment 2 because w(po, c) > — 1 
and ycosx = — y < — 1 on segment 2. Thus, if \I/ n (&) = $(c) for some 6 e (0,6^) and some c, then the functions 
u(p, c) and U(p, b) and their derivatives match at po, hence we have a solution defined on the whole interval < p < 1 
in the n th nodal class. The integer part of (6(1,6) — 6(0, c))/tt counts the number of zeros of w. In particular, 
*o(0) = $(0) since u(p, 0) = U(p, 0) = is a solution in the th nodal class. 

Remark: The choices of 0(1/2,0) = n and 0(1,0) — ir were made so that the indexing of solutions is correct. 
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FIG. 3: Illustration to the proof of Proposition 3. 

Now we are ready to prove: 

Proposition 3. For any positive integer n there is B n < boo such that the solution U(p, B n ) is in the n th nodal class. 

Proof: Since boo < bo, there is an po < 1 such that U(po, b) < bo for all b <E [0, boo]. Note that U'(l, b) < for b < b^ 
and thus U'(p, b) < on p < p < 1 because U"(p, b) > if U'(p, b) = and U(p, b) < b Q . It follows that U(p, b) > 
and hence W(p,b) > — 1 for po < p < 1 and < b < boo. From this conclude that the curve ^ n (b) is simple, i.e., 
^n(b) = ^ n (b) implies that b — b. For nonnegative b and b this follows from the uniqueness of solutions of ODEs, 
while for b < < b, , 5 rl (6) = $„(()) is impossible because W(po,b) > —1 by choice of po and ycosx — ~y < — 1 
on segment 2 of the curve. Since lim^f^ R(po,b) = -R(/9o,^oo) = and limb^-oo R(po,b) = oo, it follows from the 
Jordan curve theorem that the curve ^ n separates the half plane into two open regions: i?5_\image(4' n ) = AuB. 
Moreover, points p and q are in the different components A and B iff there is a curve from p to q crossing image(^> n ) 
transversally at exactly one point. 

Note that 9(1, b) > tt/2 if b < boo and 6'(p, b) > when 6(p, b) = tt/2, hence Q(p, b) > vr/2 for b < b^. The region 

Z = {(x,y) E R 2 + \x <ir/2-2mr} (69) 

is connected and does not meet image(^ n ) because Q(p,b) > 7r/2, hence Z must be contained in A or B\ suppose 
that Z C A. The curve $ meets A because lim^oo &(p, c) = — oo by Lemma 2 and hence the curve <& meets Z C A. 
Note that image(^ n ) D image(^> m ) = if n ^ m by uniqueness of solutions of ODEs, hence ^(O) ^ image(^> n ) if 
n > 0. 

We will now show that $(0) G B. That will complete the proof because the curve $ is connected, hence it must 
cross image^n), which means that there is a solution in the n th nodal class. To show this we construct a curve in 
R^_ from $(0) to Z that crosses imaged n ) transversally at exactly one point. We denote by M an upper bound of 
R{po, b) on < 6 < boo and define a curve 

= f (tt, 1 + t) for < t < M (segment 1), 

](n + M-t,M+l) for M < t < (2n + l)ir (segment 2). 

Note that segment 1 of T cannot meet segment 1 of ^ n because segment 1 of ^ n has w > — 1 and segment 1 of T 
has ycosx < —1. Also, segment 1 of T cannot meet segment 2 of ^ n because they have different first coordinates 
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ir — 2ir < ir. Finally, segment 2 of T cannot meet segment 1 of ^> n because R(po, b) < M < M = 1. Thus, the only 
intersection point of T and ^f n is at (tt — 2mr, M + 1). This completes the proof of Proposition 3. 

Proposition 4. For any positive integer n there is B n > such that the solution U(p, b n ) is in the n th nodal class. 

Proof: The proof proceeds along the same lines as the proof of Proposition 3 except for some technical modifications 
which we leave to the reader. 

Theorem 1 is the immediate consequence of Propositions 3 and 4 if we set 6271-1 = B n and bi n = B n . 
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